arXiv:1508.03347vl [astro-ph.HE] 13 Aug 2015 


Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 17 August 2015 (MN style file v2.2) 


Asymmetric supernova remnants generated by Galactic, massive 
runaway stars 

D. M.-A. Meyer,N. Langer,^ J. Mackey,^’^ P. F. Velazquez^ and A. Gusdorf^ 

^Argelander-Institut fUr Astronomie der Universitdt Bonn, AufdemHUgel 71, 53121, Bonn, Germany 
^1. Physikalisches Institut, Universitdt zu Kdln, ZUlpicher Strafie 77, 50937 Kdln, Germany 
^Instituto de Ciencias Nucleares, UN AM, Apartado Postal 70-543, 04510 Mexico, DF, Mexico 

"^LERMA, UMR 8112 du CNRS, Observatoire de Paris, Ecole Normale Superieure, 24 rue Lhomond, 75231 Paris Cedex 05, France 


Received January 18 2015; accepted Month day, 2015 


ABSTRACT 

After the death of a runaway massive star, its supernova shock wave interacts with the 
bow shocks produced by its defunct progenitor, and may lose energy, momentum, and its 
spherical symmetry before expanding into the local interstellar medium (ISM). We investigate 
whether the initial mass and space velocity of these progenitors can be associated with asym¬ 
metric supernova remnants. We run hydrodynamical models of supernovae exploding in the 
pre-shaped medium of moving Galactic core-collapse progenitors. We find that bow shocks 
that accumulate more than about 1.5 Mq generate asymmetric remnants. The shock wave first 
collides with these bow shocks 160 — 750 yr after the supernova, and the collision lasts until 
830 — 4900 yr. The shock wave is then located 1.35 — 5 pc from the center of the explosion, 
and it expands freely into the ISM, whereas in the opposite direction it is channelled into the 
region of undisturbed wind material. This applies to an initially 20 Mq progenitor moving 
with velocity 20 km and to our initially 40 Mq progenitor. These remnants generate mix¬ 
ing of ISM gas, stellar wind and supernova ejecta that is particularly important upstream from 
the center of the explosion. Their lightcurves are dominated by emission from optically-thin 
cooling and by X-ray emission of the shocked ISM gas. We find that these remnants are likely 
to be observed in the [Olll] A 5007 spectral line emission or in the soft energy-band of X-rays. 
Finally, we discuss our results in the context of observed Galactic supernova remnants such 
as 3C391 and the Cygnus Loop. 
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1 INTRODUCTION 


Massive stars are rare but crucial to understand the cy¬ 
cle of matt e r in the interstellar medium (ISM) of galax¬ 
ies ( Langen l2012r). Significantly influenced by their rota - 
tion ( Langer et al.lll999l: Ivan Marie etal. 1120081 : IChita et Zll2008h . 
bulk motion ( Brighenti &^^Er colel 1 1995'^ ) or by the pres¬ 
ence of a companion (IStevens et al.l Il992h . their strong winds 
shape their s urroundings and chemically augment their ambi¬ 
ent medium (IVinkI l2006h . Some of these stars explode as lu¬ 
minous supernovae whi ch release ejecta interacting with their 
pre-shaped environmen t ( Borkowski et all 1 19921 : IVink et all 1 19961 . 


I 1997 I : l^n Veelen et al.ll2()09l) . This event gives birth to supernova 

remnants replenishing the ISM with momentum and kinetic energy 
up to about 120 pc from the center of the explosion (iBadenes et al.l 

l2^ . 


* E-mail: dmeyer@astro.uni-bonn.de 


Supernovae have been noticed in ancient Asia with the 
naked eye, e .g. the guest-star recorded in AD 185 by Chinese 
astronomers (iGreen & Stephensonl l2003h . whereas the first su¬ 
pernova remnant ha s been ident ified spectroscopically almost 
two millennia later (iBaadd Il938h . Nowadays, surveys provide 
us with obs ervations of Gal actic su pernova remnants, e.g. in 
gamm a-rays ( Abdoet^ 201Q|) . X-ray JPannuti et al.l2014 ), the in¬ 
frared (iReach et al.ll2006 ) or in submillimeter (Ivan Dishoeck et aP 
I 1993 I) . Catalogues of remnants visible in the radio waveband in 
the no rther n and southern hemisphere are available in iKothes et al.l 
(l2006h and lWhiteoak & GreenI d 19961) . respectively. An exhaustive 
ca talogue of the known Galactic supernova remnants was compiled 
bylo reenl (l2009l) . These abundant observations reveal a diversity of 
complex morphologi es such as shell s, annuli, cylinders, rings and 
bipolar structures (cf. lGaensle^l 19991) . 


The shape of young supernova remnants depends (i) on the 
geometry of the supernova ex plosion and (ii) on the (an)isotropy 
of their ambient medium (IVinkI l2012l) . They therefore ex- 
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hibit a wide range of morphologies that can be used to con¬ 
strain their progenitors and/or ambient medium properties. The 


properties 


borders of neighbouring diffuse nebulae or filaments that affect the 
propagation of the supernova ejecta. Models o f remnants develop¬ 
ing in a pre-e xisting wind cavity are shown in lTenorio-Tade et~^ 
(1 19901 Il99lh . and demonstrate that mixing of material happens 
in the former wind bubble. Multi-dimensional models of the for¬ 
mation o f knots by wind-wind collis ion around Cassiopeia A are 
shown in IPerez-Rendon et all (l2009h and the effects of this frag¬ 
mented Wolf- Rayet shell on the rebri ghtening of young remnants 
is explored in Ivan Veelen et al.l (l2009l) . Supernova remnants devel¬ 
oping through an edge of a dense region, e.g. a molec ular cloud, 
give rise to champagne fiows jTenorio Tade et alll985h . If the su¬ 
pernova happens near a denser region, the reverse shock is reflected 
towards the center of the explosion and a hot region of shocked ma¬ 
terial forms (iFerreira & de Jage3l2008h . A strong magnetization of 
the ISM can induce a co llimation of the supernova ejecta e ngender- 
ing elongated remnants (iRozvczka & Tenorio-Tadeiri995h . 

Particularly, the bow shocks produced by runaway massive 
stars are an ideal site for the generation of an anisotropic cir- 
cumstellar distribution. This is likely to happen in the Galactic 
plane, where most of the massive stars w hich are bo t h in the 
field and classified as runaw ay are found (iGie 3 1 1981 E aauwl 
1 19931 : iHuthoff & Kapej l2002h . A few of them are identified as 
evolved massive stars and three of them ar e red supergiants with 
detec te d bow sho c ks, i. e. Betelgeuse ( Noriega-fiespo et"aD 


_ ga _ 

1 19971 : IPecin et ^ l2012l), ii Cep dCox et al.l 2012 ) and 


IRC—10414 ( Gvaramadze et al.l l2014l) . Consequently, and 
because these stars will explode as core-collapse supernova, their 
circumstellar medium is of prime interest in the study of aspherical 
supernova remnants. 

The circumstellar mediu m of Galactic runaway red super giant 
stars is numerically studied in iBrighenti & D’Erc^ (1 1995 al io) as 
an attempt to explain non- sphe rical supernova rem nants. The works 
by Ivan Marie et al.l (|201 ll) and lPecin et al.l (l2012h tailor models to 
Betelge use’s bow shock and estimate in the context of recent obser¬ 
vations jCox et al.l2012h how the drag force on dust grains modifies 
the evolution of its contact discontinuity. The effects of the mass 
loss and space velocity on the shape and lumin osity of bow shocks 
around red supergiant stars is investigated in iMever et aP (l20l4 
hereafter Paper I). The repercussions of a weak ISM magnetic field 
on the damp ing of instabilitie s in th e bow shocks of Betelgeuse is 
presented in Ivan Marie et al.l (l2014l) . The stabilizing role of pho¬ 
toionization by an external source of rad i ation on the bow shock 
of IRC—10414 is shown in iMever et al.l (l2014h . The above cited 
models can be understood as investigations of the circumstellar 
medium of Galactic runaway core-collapse progenitors near their 
pre-supernova phase. 

After the supernova explosion, the forward shock 
of the blastwave interacts w i th the free-stre aming stellar 
wind jChevalier & Liangl Il989l : IChevaliej Il982l) . Supernovae 
showing evidence of interaction with circumstellar structures are 
commonly denoted as Type Iln and their corresponding lightcurves 
provide infor mation on the progenitor and properties of its close 
surroundings jSchlegell I199QI : iFilippenk^ 1 19971 : Ivan Marie et al.l 
bold) . About 10—lOOyr after the explosion, the shock wave 
collides with the bow shock along the direction of motion of its 
progenitor, wherea s it expands in a cavity of wind material in the 
opposite direction (iBorkowski et al.ll 19921) . 


In this spirit jRozvczka et d] (1 19931) model supernovae in oval 
bubbles generated by moving progenitors. They neglect the progen¬ 
itor stellar evolution but demonstrate that elongated jet-like struc¬ 
tures of size of about 10 pc form when the shock wave expands 
into the wind bubble. A model interpreti ng the cool iet-like [O llll 
A 5007 feature found in the Crab nebula JBlandford et as 

a shock wave channelle d into the t r ail pro du ced by its progenitor’s 
motion is presented in ICox et akl (Il99ll) . IBrighenti & D’Erc^ 
(1 19941) show that if the runaway progenitor evolves beyond the 
main-sequence phase, the supernova explosion happens out of the 
main-sequence wind bubble, and the subsequent remnant develops 
as an outfiow upstream from the direction of motion of the progen¬ 
itor. 

In this work, we aim to determine the degree of anisotropy of 
supernova remnants generated by runaway core-collapse progeni¬ 
tors moving through the Galactic plane. We model the circumstel¬ 
lar medium from near the pre-supernova phase of a representative 
samp le of the most common runaway massive stars (lEldridge et al.l 
I2OIII) . We calculate one-dimensional hydrodynamic models of the 
supernova ejecta interacting with the stellar wind and use them as 
initial conditions for two-dimensional simulations of the supernova 
remnants. Finally, we discuss the emitting properties of the most 
aspherical of these remnants. 

This project is diff erent from previous stud¬ 
ies (IBrighenti & D’ Ercol^l 1 994l) because (i) we use self-consistent 
stellar evolution models, (ii) we consider both optically-thin 
cooling and heating along with thermal conduction, (iii) we trace 
the mixing between ISM, stellar wind and supernovae ejecta 
inside the remnants and (iv) our grid of models explores a broader 
space of parameters than works tailored to a particular supernova 
remnant, e.g. Kepler’s s upernova remna n t dBorkows ki et al.ll 19921: 
Velazquez ^ al.l I2OO6I : IChiotellis et al.l I2OI2I: iToledo-Rov et^ 
2 OI 4 I) or Tycho s supernova remnant dVigh et al]l201 ll) . We neglect 
the magnetization, inhomogenity and turbulence of the ISM and 
ignore the cooling in the sh ock wave induced by the production 
of G alactic cosmic rays dOrlando et al.l I2OI2I : ISchure & Belli 
I2OI3I) . We assume that the supernova explosions do not have 
any intrinsic anisotropy. Furthermore, we assume than no pulsar 
wind nebula remains inside the supernova remnants and modifies 
the refiect ion of the reverse sh ock towards the center of the 
explosion dBucciantini et al.ll2003h . 

This paper is structured as follows. We begin Section[2|by dis¬ 
cussing our numerical methods and initial parameters. The mod¬ 
elling of the circumstellar medium of our progenitors is shown in 
Section 13 We describe the calculations of supernova remnants de¬ 
veloping inside and beyond their progenitors’ bow shock in Sec¬ 
tions |4] and |5] respectively. Section [6] discusses and compares our 
models of aspherical remnants with observations. We conclude in 
Section |7] 


2 METHOD AND INITIAL PARAMETERS 

In this section, we review the numerical methods used to model 
the circumstellar medium of our progenitors and we present the 
procedure to set up supernova blastwaves. 


2.1 Modelling the circumstellar medium 

We per fom two-dimensiona l sim ulations using the code 
PLUTO dMignone et al.l 120071 I 2 OI 2 I) to model the circumstel¬ 
lar medium of moving core-collapse supernova progenitors. We 
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Asymmetric, Galactic supernova remnants 3 


Table 1. Wind properties at the end of the used stellar evolution mod¬ 
els, at fpsn. Mi, (Mq) is the initial mass of each star, M (M 0 yr“^) 
their mass loss and u^ind (kms“^) their wind velocity at a distance of 
0.01 pc from the star, resp ectively. Teff (K) is the effective temperature of 
the stars iBrott et al.ll201 ih . 


gas, where r is the position vector of a grid cell. Its value is set to 
Qi{r) = 1 for the wind material and to Qi(r) = 0 for the ISM 
gas. Our cooling curves and numerical methods are extensively de¬ 
tailed in section 2 of Paper L Results are presented in Section [3] 


M„ (Mq) 

M(Moyr-l) 

rwind(kms 1) 

Teff (K) 

20 

10 - 6.11 

16 

3.2 X 10^ 

40 

10-4.79 

11 

3.2 X 10^ 

70 

10-4.48 

50 

5.8 X 10^ 


solve the equations of hydrodynamics in a cylindrical computa¬ 
tional domain (O; R, z) of origin O, which is coincident with the 
location of the runaway star and has rotational symmetry about 
R = 0. A uniform grid of Wr x Nz grid cells is mapped onto a 
domain of size [0, i^max] x [zmin, ^max], respectively. We define R 
and z as the unit vectors of the axis OR and Oz, respectively. The 
gri d spatial resolution is A = i?max/WR. Following the method 
of IComeron & Kap^ (Il998h . we release the stellar wind on a 
circle of radius 20 grid cells centered on the origin and compute 
the wind-ISM interaction in the frame of reference of the moving 
progenitor. 

We model the circumstellar medium of initially 10, 20 and 
40 M© stars moving with space velocity ranging from u* = 20 
to 70kms“^. The (time-dependent) stellar wind pro perties are 
taken from stellar evolution models terott et al.ll20TIh . We con¬ 
sider a hom ogeneous ISM with a hydrogen number density nu = 
0.57 cm“^ dWolfire et al.ll2003h . i.e. we assume that the stars are 
exiled from their parent cluster and move in the low-density ISM. 
We set the ISM temperature to Tism ~ 8000 K and include 
gain/losses by optically-thin radiative cooling assuming that the gas 
has solar metallicity (sections 2.3 and 2.4 of Paper I ) All our bow 
shock models include e lectronic thermal conduction dspitzeill 19621 : 
ICowie & Mc^ll977h . 

We start our models at tstart ~ tpsn — 32 tcross, where tpsn is 
the time at the end of the stellar evolution model. The time inter¬ 
val 32 tcross is sufficient to simulate these bow shocks getting rid 
of any switch-on effects arising during their development, where 
Across = R{0)/vi, is the cr ossing-time and R(0) is the stand¬ 
off di stance of the bow shock (iBaranov. Krasnobaev & Kulikov ski j 
Il97lh . The calculation of each bow shock model is followed until 
the end of the stellar evolution model, at tend- Note that our ini¬ 
tially 10 and 20 M© progenitors explode as a red supergiant (Pa¬ 
per I). We run these models with the same underlying assumptions 
as in Paper I, especially consideri ng that the stefiar ra diation field 
is not trapped by the bow shock dWeaver et al.ll 19770 . Moreover, 
we assume that the evolutionary model of our initially 40 M© pro¬ 
genitor, that does not go all the way up to the pre-supernova phase, 
is sufficient to approximate the mass distribution at the time of the 
supernova explosion. Our stellar evolution models are described in 
section 2.2 of Paper I and the wind properties at tpsn are shown in 
TablelU 

Our models of the circumstellar medium from near the pre¬ 
supernova phase are named with the prefix PSN followed by the 
initial mass of the progenitor M* (first two digits, in M©) and its 
space velocity u* (two last digits, in kms“^). We adopt grid di¬ 
mensions such that it includes the wake of the bow shocks produced 
during about u* lOtcross (our Table[2l) in order to properly model the 
expansion of the shock wave through the tail of the bow shock up to 
times of the order of 10^ yr. The stellar wind is distinguished from 
the ISM material using a scalar Qi (r) passively advected with the 


2.2 Setting up the supernova shock wave 


We perform one-dimensional hydrodynamical simulations of the 
shock wave expanding into the stellar wind. The blastwave is char¬ 
acterized by its energy fixed to Eei = 10^^ erg and by the mass of 
the ejecta Mej. The latter is estimated as, 

r'tpsn 

Mej =Mi,- M{t)dt - Mco, (1) 

Jto 

where to and tpsn are the time at the beginning and the end of 
the used stellar evolution model, respectively. Note that we assume 
Mej for our 40 M© progenitor, because we ignore its post-main- 
sequence evolution. The quantity Mco = 2 M© in Eq. 0 is the 
assumed mass of the residual compact object left after the super¬ 
nova (our Table [3]). 

We set up the super nova using the meth od detailed 
in I Whalen et al.l (l2008l) and in Ivan Veelen et al.l (l2009h . It assumes 
that the blastwave density profile p{r) is a radial piece-wise func¬ 
tion of the distance r to the center of the explosion in the region 
[0, Tmax], where rmax is the radius of the shock wave when we start 
the simulations. Under these assumption, the ejecta density profile 
is, 

{ Pcore (j ') 

Pmax(r’) 

Pcsm {t ) 

where, 

1 (ios7^)-®/y_3 

Pcorelrj - (3^n-3)_5/2 

is constant up to the inner core of radius rcore and, 

1 / r 

Pmax(r) - (3JV^n-3)(n-5)/2 ( 

is a steeply decreasing function of inn er radius rcore and external 
radius rmax (iTruelove & McKe^l 19991) . The power law index n of 
Eq. (H-dS is set to the us ual value n — 11 for core-collapse 
supernovae (IChevalieilll982l) . In relation (O, pcms is the freely- 
expanding wind profile measured from the simulations along the 
symmetry axis Oz, in the direction of motion of the progenitor 
{z ^ 0). We use it as initial condition in the [rmax,'rsncsm] of the 
domain, where rgncsm < 77(0) is outer border of the domain. 

The ejecta obey a homologous expansion, i.e. the velocity pro¬ 
file v(r) is. 



if r ^ r core 5 

if r core r r” max 5 (2) 

if r ^ r max 5 


r 

r(r) = if t > tmax, (5) 

where t is the time after the supernova explosion. The ejecta veloc¬ 
ity at rcore is therefore. 


/ 10(n-5)E;ej \ 


(ITruelove & McKe^l 19991) . The choice of rmax is free, as long as a 
mass of stellar wind smaller than Mej is enclosed in [rmax, rgncsm]- 
W e determined i ts rmax using the numerical procedure described 
in I Whalen et al.l (l2008l) . We start the simulation at tmax = 


© 0000 RAS, MNRAS 000, 000-000 















































4 D. M.-A. Meyer et al. 


Table 2. Input parameters used in our simulations of the bow shocks generated by supernova progenitors. Parameters and i^max are the limits of the 
domain along the a:-axis and y-axis (in pc), respectively. and Nz are the number of cells discretising the corresponding directions and A is the grid 
resolution (in pc cell“^). The simulations start at a time fstart (in Myr) after the star’s birth and are run until the end of the used stellar evolution models, at 
fpsn (in Myr). 


Model 

M* (Mq) 

i;*(kms ^) 

-2^111111 (pc) 

-Rmax (pc) 

Nk 

Nz 

Z\ (pc cell ^) 

f start (Myr) 

fpsn (Myr) 

PSN1020 

10 

20 

-6.0 

7.0 

1225 

1400 

5.71 X 10-3 

23.7 

24.7 

PSN1040 

10 

40 

-2.1 

2.0 

1143 

1600 

1.75 X 10-3 

24.5 

24.7 

PSN1070 

10 

70 

-0.9 

0.8 

950 

1425 

8.42 X 10-"^ 

24.6 

24.7 

PSN2020 

20 

20 

-30.0 

35.0 

2333 

2333 

1.49 X 10-2 

6.63 

9.05 

PSN2040 

20 

40 

-9.0 

8.0 

600 

1200 

1.33 X 10-2 

8.65 

9.05 

PSN2070 

20 

70 

-4.5 

4.0 

1067 

1600 

3.75 X 10-3 

8.70 

9.05 

PSN4020 

40 

20 

-90.0 

100.0 

1500 

1575 

6.00 X 10-2 

0.0 

4.5 

PSN4040 

40 

40 

-36.0 

30.0 

1500 

2100 

2.00 X 10-2 

2.00 

4.5 

PSN4070 

40 

70 

-30.0 

20.0 

1000 

2000 

2.00 X 10-2 

3.50 

4.5 


Table 3. Simulations parameters used in our simulations of supernovae in¬ 
teracting with the unperturbed stellar wind. Parameter Mej is the mass of 
the ejecta (in Mq) and rmax is the size of the one-dimensional spherically 
symmetric domain (in pc). The simulations are started at f = 0.04 yr. The 
last column indicates the time at the end of our simulations, fsncsm (in yr). 


Model 

Mej (Mq) 

rmax (pc) 

fsncsm (yr) 

SNCSM 1020 

7.7 

0.30 

40 

SNCSM 1040 

7.7 

0.20 

13.5 

SNCSM 1070 

7.7 

0.13 

15 

SNCSM 2020 

17.7 

0.90 

40 

SNCSM 2040 

17.7 

0.50 

25 

SNCSM 2070 

17.7 

0.25 

20 

SNCSM 4020 

20 

3.00 

400 

SNCSM 4040 

20 

1.50 

200 

SNCSM 4070 

20 

0.90 

180 


rmax/ t;max, whcrc t;max is sct to 30000 kms ^ (Ivan Veelen et al.l 

l2009l> . 

We choose a uniform grid of resolution A ^ 10 ^ pc cell ^ 
and follow the expansion of the shock wave until slightly before it 
reaches the reverse shock of the bow shocks produced by our pro¬ 
genitors. These models are labelled with the prefix SNCSM (our 
Table [3]). Additionally, we use a second passive scalar Q 2 {r) to 
distinguish the ejecta from the stellar wind. We carry out these 
one-dimensional calculations using a uniform spherically symmet¬ 
ric grid. We use a finite volume method with the Harten-Lax-van 
Leer approximate Riemann solver, and integrate the Euler equa¬ 
tions with a second order, unsplit, time-marching algorithm. Dis¬ 
sipative processes are computed using our cooling curve for fully 
ionized gas. Results are presented in Section lU 


2.3 Modelling the supernova remnants 

In order to resolve both the early interaction between the blastwave 
interacting with the circumstellar medium and the old supernova 
remnant, we adopt a mapping strategy. We run two-dimensional hy- 
drodynamical simulations of the shock waves interacting with the 
bow shocks using a squared computational domain of size about 
4i?(0) which is supplied with a uniform rectilinear grid. These 
models are labelled with the prefix YSNR. The above-described 
one-dimensional simulations of the ejecta interacting with the stel¬ 


lar wind are mapped into a circle of radius rmax < ^(0) centered 
on the origin O of the domain. We run these simulations starting at 
tsncsm until thc shock wave has passed through the forward shock 
of the bow shock and reaches a distance of about 2i?(0) in the di¬ 
rection of motion of the progenitor, at tysnr (our Table ID). 

The remnants at tysnr are mapped a second time onto a 
larger computational domain which includes both the entire pre¬ 
calculated circumstellar medium and the calculations of the young 
supernova remnants (our Tables [2| and 01 ). The regions of this do¬ 
main which overlap neither the bow shock nor the remnant are filled 
with unperturbed ISM gas. We start the simulations at time tysnr 
and follow edge of the domain in the —z direction, at tosnr. These 
simulations are labelled with the prefix OSNR (our Table [S]). Re¬ 
sults are presented in Section [5] 


3 THE PRE-SUPERNOVA PHASE 


In Fig. [T] we show the gas density fields in our bow shocks from near 
the pre-supernova phase in the models PSN1020 (initially 10 Mq 
star, Vi, — 20kms“^, Fig. [1^), PSN1040 (initially 10 M© star, 
Vi, = 40kms“^, Fig. [H)) and PSN1070 (initially 10 M© star, 
Vi, = 70kms“^, Fig. [T]:). Figs. [2] and [3] show the same, but for 
our 20 and 40 M© stars. The figures plot the density at a time tpsn 
and do not show all of the computational domain. In Figs.[T]|2]and[3] 
the overplotted solid black line is the contact discontinuity, i.e. the 
border between the wind and ISM gas where the value of the pas¬ 
sive scalar Qi{r) = 1/2. 

The bow shocks of our 10 and 20 M© stars have mor¬ 
phologies consistent w ith previous studies dvan Marie et al JImll; 
iMohamed et al.l 12012, Paper I). Th eir overall shape is rather sta¬ 
ble if Vi, ^ Uw ( Dgani et^ 1 19961) and the flow inside the bow 
shocks is laminar (Fig. [T^). The bow shocks are unstable and 
exhibit Rayleigh-Taylor and/or Kelvin-Helmholtz instabilities for 
Vi, ^ Uw because of the large density difference between the 
dense red supergiant wind and the ISM gas (Figs. [I]3,c andia,b). 
For high-mass stars moving with large space velocities, e.g. the 
models with Mi, ^ 20 M© and u* ^ 40kms ~^, the shocked 
layers develop non-linear thin-shell instabilities ( Vishniad 1 1994 


iGarcia- Segura et al.]| 19961 : iBlondin & Koerwedl 19981) and induce j 

strong mixing in the wakes of the bow shocks (Fig.|2t). 

The stellar motion displaces the position of the star 
fro m the center of the cavi ty of unshocked wind mate¬ 
rial dBrighenti & D’Ercolel 1995 a lia), and this displacement is larger 
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Table 4. Input parameters used in our simulations of the supernova blastwaves interacting with the bow shocks of our progenitors. As input we use the solution 
of the shock waves interacting with the stellar winds (our Table [3. The grid parameters are similar as in our Table |2] Our simulations start at fsncsm (our 
Table|2) shortly before that the shock wave interacts with the bow shock and the models are run until fysnr (in yr) once the shock wave has gone through it. 


Model 

Input 

^^min (pc) 

Rmayi (pc) 

Nk 

iVz 

Z\ (pc cell ^) 

fsncsm (yr) 

fysnr (yr) 

YSNR1020 

SNCSM1020 

-2.0 

2.0 

1000 

2000 

2.0 X 10-3 

40 

264 

YSNR1040 

SNCSM1040 

-1.3 

1.3 

1000 

2000 

1.3 X 10-3 

13.5 

150 

YSNR1070 

SNCSM1070 

-0.7 

0.7 

1000 

2000 

4.0 X 10-4 

15 

60 

YSNR2020 

SNCSM2020 

-8.0 

8.0 

1000 

2000 

8.0 X 10-3 

40 

2400 

YSNR2040 

SNCSM2040 

-2.0 

2.0 

1000 

2000 

2.0 X 10-3 

20 

450 

YSNR2070 

SNCSM2070 

-1.0 

1.0 

1000 

2000 

1.0 X 10-3 

25 

170 

YSNR4020 

SNCSM4020 

-25 

25 

1000 

2000 

2.5 X 10-2 

400 

4900 

YSNR4040 

SNCSM4040 

-7.0 

7.0 

1500 

3000 

4.7 X 10-3 

200 

1360 

YSNR4070 

SNCSM4070 

-4.0 

4.0 

1500 

3000 

2.7 X 10-3 

180 

830 


Table 5. Input parameters used in our simulations of the supernova blastwaves interacting with the tails of the bow shocks generated by our progenitors. As 
input we use the solution of the shock waves interacting with the bow shocks (our Table |4)- The grid parameters are similar as in our Table [2| Our simulations 
start at fysnr (our Table|4) and the models are run until fosnr (in yr). 


Model 

Input 

^min (pc) 

-Rmax (pc) 

Nk 

Vz 

A (pc cell ^) 

fysnr (yr) 

fosnr (yr) 

OSNR1020 

YSNR1020 

-6.0 

6.0 

500 

1000 

1.2 X 10-2 

264 

1500 

OSNR1040 

YSNR1040 

-2.1 

2.1 

500 

1000 

4.2 X 10-3 

150 

1300 

OSNR1070 

YSNR1070 

-0.9 

0.9 

500 

1000 

4.0 X 10-"^ 

60 

1300 

OSNR2020 

YSNR2020 

-30.0 

25.0 

1000 

2000 

9.0 X 10-2 

2400 

21100 

OSNR2040 

YSNR2040 

-9.0 

9.0 

1000 

2000 

2.5 X 10-3 

450 

15000 

OSNR2070 

YSNR2070 

-4.5 

4.5 

1000 

2000 

9.0 X 10-3 

170 

10000 

OSNR4020 

YSNR4020 

-90.0 

70. 

1000 

1714 

7.0 X 10-2 

4900 

49500 

OSNR4040 

YSNR4040 

-35.0 

25.0 

1000 

2200 

2.5 X 10-2 

1360 

14000 

OSNR4070 

YSNR4070 

-30.0 

15.0 

700 

2333 

2.1 X 10-2 

830 

10500 


for velocities t;* ^ 20kms“^ (Fig. |2t). The bow shocks which 
have the most pronounced tunnels of low-density gas are produced 
either by our initially 20 M© star moving with 20 km s“ ^ or by our 
initially 40 M© star (Fig.|2t and[^-c). In the region downstream 
from the progenitor, the reverse shock, which forms the walls of 
the cavity, has a rather smooth appearance for u* ^ 20kms“^ 
(Fig. [2^) but it is ragged for u* ^ 40kms“^ (Fig. [3:). Finally, 
note that the model PSN2020 has a double bow shock due to the fi¬ 
nal increase of the mass loss that ends the red supergiant phase. 
This structure is called a Napoleon’s hat and it develops when 
the bow shock from a new mass-loss event go es through the one 
generated by the p revious evolutionary phase (IWang et al.lll993l : 
iMackev et al.ll201 A . 

The stand-off distance R{0) and the mass M trapped in the 
bow shocks upstream from the star (z ^ 0 ) are summarised in 
Table [6| The more massive bow shocks are the biggest ones, e.g. 
our bow shock model PSN4020 has the largest stand-off distance 
R{0) 5pc and has accumulated about 116 M© of shocked gas. 
They are generated by high-mass stars moving with small space 
velocities, i.e. M* ^ 20 M© and u* ^ 40kms“^ (Figs.|2^-|3^b). 
In Fig. 0] we show the average density profiles in our simulations 
of our 10 (Fig.|4^), 20 (Fig.ldJ)) and 40 M© models (Fig.|4]:), that 
we use as initial conditions for our one-dimensional simulations of 
the supernova shock waves interacting with their surroundings, see 
Eqs. (ID)—13). Fig. 01 illustrates that the most massive bow shocks 
have the largest R{0), i.e. they are the most voluminous and are 


Table 6. Gas mass M (M©) in the region of the bow shocks that is up¬ 
stream from the progenitor (z ^ 0) and stand-off distance at the contact dis¬ 
continuity (in pc) measured in our simulations of the circumstellar medium 
from near the pre-supernova phase, at a time f psn • 


Model 

M(M0) 

R{0) (pc) 

PSN1020 

0.06 

0.42 

PSN1040 

0.03 

0.25 

PSN1070 

0.01 

0.17 

PSN2020 

3.87 

1.35 

PSN2040 

1.10 

0.64 

PSN2070 

0.75 

0.55 

PSN4020 

116.00 

5.00 

PSN4040 

9.40 

2.70 

PSN4070 

1.65 

1.55 


reached by the shock wave about R(0)/v{r) ^ lO^yr after the 
explosion (Section01). 


4 THE YOUNG SUPERNOVA REMNANT 

This section presents our simulations of the supernova blastwaves 
interacting with their circumstellar medium until the shock wave 
reaches the outer layer of its surrounding bow shock. Spherical 
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Figure 1. Grid of stellar wind bow shocks from the pre-supernova phase 
of our initially 10 M© progenitor as a function of its space velocity with 
respect to the ISM, with velocity u* = 20 (a), 40 (b) and 70kms“^ (c). 
The nomenclature of the models follows Table[2| The bow shocks are shown 
at fpsn. The gas number density is shown with a density range from 10“ ^ to 
30.0 cm“^ in the logarithmic scale. The white cross marks the position of 
the star. The solid black contour traces the boundary between wind and ISM 
material Qi{r) = 1/2. The x-axis represents the radial direction and the 
y-axis the direction of stellar motion (in pc). Only part of the computational 
domain is shown in the figures. 

remnants are distinguished from asymmetric remnants as a func¬ 
tion of their progenitors’ properties. 

4.1 The ejecta-stellar-wind interaction 

In Fig. [5] we show a typical interaction between a supernova shock 
wave and the surrounding stellar wind, before the shock wave col¬ 
lides with the bow shock. We plot the gas number density (solid 
blue line) and temperature (dashed red line) profiles in our model 
SNCSM1020 at a time about 22 yr after the supernova explosion. 


Figure 2. As Fig. [T] with our initially 20 Mq progenitor. Note the 
Napoleon’s hat structure of the bow shock in panel (a). 

It assumes a release of Mej = 7.7 M© of ejecta together with a 
kinetic energy of Eei = erg (our Table |3). The structure is 
composed of 4 distinct regions: the expanding ejecta profile, itself 
made of two regions, the core and the envelope (iTruelove & McKe3 
Il999h . the shell of swept up shocked ejecta and shocked wind ma¬ 
terial , and finally the undisturbed circumstellar material (iChevalied 
Il982h . 

The shell is bordered by two shocks, a reverse shock that is 
the interface between ejecta and shocked ejecta, and a forward 
shock constituting the border betwee n shocked wind and undis- 
turbed freely-expanding stellar wind (IChiotellis et al.ll201^ . The 
core of the ejecta (r < 0.08 pc) has a v ery low temperature 
because its thermal pr essure is initially null (IWhalen et al.ll2008l : 
Ivan Veelen et al.ll200^ . The temperature slightly increases up to a 
few tens of degrees in the envelope of ejecta (0.08 ^ r ^ 0.16 pc) 
because (i) we use a homologous velocity profile which results in 
increasing the thermal pressure close to the high-velocity shock 
wave and (ii) the decreasing density pmax oc increases the 
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Figure 3. As Fig.[T] with our initially 40 Mq progenitor. 



temperature T oc p/p. At radii r 0.16 pc is the hot (T ^ 10^ K) 
and dense (n ^ 10^ cm“^) gas. This region between the shell and 
the shock wave is hot because it is shock-heated by the blastwave 
and it has not yet cooled. All our models have a similar behaviour. 




Figure 4. Gas number density profiles (in cm taken along the direction 
of motion of our initially 10 (a), 20 (b) and 40 M© (c) progenitors at fpsn- 


4.2 The shock wave interacting with the bow shock 

The supernova interacting with the bow shock generated by the 
10 M© star moving with u* = 40kms“^ is illustrated in Fig.[6l 
We show the density field in our simulation YSNR1040, composed 
of a shock wave interacting with its circumstellar medium (our 
model SNCSM1040). The density stratification is shown at times 
tsncsm ~ 152 yr (Fig.[6^), 154, 162, 168, 192 and tysnr ~ 246 yr 
(Fig. [6f) after the explosion, respectively. Note that the part of the 
computational domain plotted in Fig.[6f is larger than in Figs. [6^- 
e. Corresponding cross-sections measured in these density fields 
along the Oz axis are shown in Fig.[71 

At a time tsncsm ~ 152 yr the solution SNCSM1040 is 


mapped on the model PSN1040 shortly before the supernova rem¬ 
nant shock wave reaches the reverse shock of the bow shock 
(Fig. it). The organisati on of the remnant is sim ilar to the situa¬ 
tion depicted in fig. 5 of I Velazquez et al.l (l2006h . From the origin 
to larger radius, the curve (a) of Fig. [7] plots the initial plateau of 
density pcore(^) ~ 10^cm“^, the steep profile pmax(r) (x r~^, 
the shell of swept-up ejecta at about 0.16 pc (cf. Fig.|5]), the shock 
wave progressing in the freely-expanding wind, the red-supergiant 
bow shock from about 0.21 to about 0.27 pc, and the unperturbed 
ISM. 

In Fig. [6j) the shock wave collides with the reverse shock of 
the bow shock and begins to interact with the dense (n 25cm“^) 
shocked wind. The interaction starts at the stagnation point be¬ 
cause this is the part of the bow shock with the smallest ra- 


© 0000 RAS, MNRAS 000, 000-000 

















































8 D. M.-A. Meyer et al 



1e+12 


- 1e+10 


- 1e+08 


- 1e+06 


10000 


- 100 


CD 

c5 

0 

Q. 

E 

<D 


Figure 5. Total number density (solid blue line, in cm“^) and temperature (dashed red line, in km s“^) profiles of the supernova ejecta interacting with the 
circumstellar medium of our initially 10 M© progenitor. Profiles are shown at about 22 yr after the supernova ex plosion. The distinct region s composing the 
supernova remnant are indicated on the figure. The substructure of the shocked ejecta follows the nomenclature in iTruelove & McKea ( Il999h . 


dius teorkowski et al.l[T992h . The blastwave decelerates and loses 
its spherical symmetry, the shock wave penetrates the reverse shock 
of the bow shock at a time 154 yr with velocity v ~ 6700 km s~^ 
whereas it hits its forward shock Syr later with velocity v ^ 
4500 km s“^. At 162 yr, the shell of shocked ejecta merges with 
the former post-shock region at the reverse shock of the bow shock, 
and its material is compressed to n ^ 85cm“^ (curve (c) of 
Fig.[7]). 

In Fig.[6l: the shock wave has totally penetrated the bow shock, 
both a reflected and a transmitted shock wave form at both the ends 
of the bow shock. In Fig.[6li the transmitted shock at the former for¬ 
ward shock starts expanding i nto the undisturbed ISM. As sketched 
in fig. 3 of iBorkowski et al.l (Il992h . a bump emerges beyond the 
bow shock and it reaches about 0.32 pc at a time 168 yr after the 
explosion (curve (d) of Fig. □i.it expands and enlarges laterally as 
the shock wave, that is no longer restrained by the material of the 
bow shock, penetrates the undisturbed ISM, accelerates and pro- 
gr essively recovers its spherica l symmetry (Fig.[6t, see also Fig. 1 
of iBrighenti & D’Ercol jl994h . 

At a time 246 yr, the shock wave has recovered its spheric¬ 
ity (Fig. [6f). Note that the shock wave is slightly constricted in 
the cavity of unshocked wind as it expands downstream from the 
direction of motion of the progenitor. This anisotropy is a func¬ 
tion of the circumstellar density distribution at the pre-supernova 
phase and governs the long term evolution of the supernova rem¬ 
nant (Section |43]) . The curve (f) of Fig. [7] shows the density struc¬ 
ture composed of the plateau whose density has diminished to 
Pcore ~ 600 cm the twice shocked ejecta, the twice shocked 
stellar wind, twice sh ocked ISM, shocked I SM and finally the 
unperturbed ISM (c.f. IBorkowski et ^Il992h . These regions are 
not clearly discernible because of the mixing at work provoked 
by the multiple reflections and refractions proliferating through 
the remnant (curve (f) of Fig. [T]). We underline that our one¬ 
dimensional, spherically-symmetric calculations do not allow us 
to model the Rayleigh-Taylor instabilities triggered at the inter¬ 
face b etween supernova eje cta and uniformly-expanding stellar 
wind jChevalier et al.l Il992h . These instabilities may affect the 
propagation of the shock wave through the bow shock. 


4.3 Spherical and aspherical supernova remnants 

In Fig. [8] we represent the density profiles of our models of super¬ 
nova remnants at a time tysnr, taken along the direction of motion 
of the progenitors. The models produced by our 10 M© progeni¬ 
tor conserve their symmetry after the shock waves collide with the 
bow shocks, e.g. our model YSNR1020 has a plateau of density 
n 10^ cm“^ at |z| ^ 1.7 pc, whereas its density distribution 
in both directions beyond the shock wave is about the ISM am¬ 
bient medium density. Here, the shock wave expansion is barely 
disturbed by the light bow shocks (our Table [6]) and the remnants 
grow quasi-spherically (Fig. [U. 

The remnant generated by our 20 M© star moving with u* = 
20 km is aspherical, in that it has a cavity of n 10“^ cm“^ 
at z < —8 pc (red curve in Fig. [SJd). This is much less pro¬ 
nounced for models with a larger u*, i.e. our models YSNR2040 
and YSNR2070, see blue and yellow curves in Fig.[8]3. They have a 
rather spherical density distribution, and only a bulge of swept-up 
wind material in the direction +£ distinguishes them from a spher¬ 
ically symmetric structure. A small accumulation of swept-up gas 
that could slightly decelerate the shock wave forms in the direction 
—z, i.e., these models do form a pronounced wind cavity (Fig.[8]3c). 

The supernova remnants of our 40 M© progenitor are all 
strongly anisotropic, e.g. our model with velocity u* = 70 kms“^ 
has a dense shell of density n 10^ cm“^ along the direction +z 
and a cavity of density n ^ 10~^ cm“^ in the opposite direction 
(Fig.[8j:). On the basis of Table[6]and according to the above discus¬ 
sion, we And that the bow shocks of runaway stars that we simulate 
and which accumulate more than about 1.5 M© generate asymmet¬ 
ric supernova remnants. We measure from our simulations that the 
collision between the shock waves and these dense bow shocks, 
located at 7^(0) 1.35 — 5 pc from the center of the explosion, 

begins about 160 — 750 yr and ends about 830 — 4900 yr after the 
supernova, respectively. In the next section, we continue our study 
focusing on the asymmetric models only, i.e. generated either by 
our 20 M© progenitor moving with u* = 20 km or produced 
by our 40 M© star. 
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Figure 6. Time sequence of a supernova interacting with the bow shock generated by our initially 10 Mq progenitor moving with velocity 40 km The 
figures correspond to times fsncsm (a) up to about fysnr (f). The gas number density is shown with a density range from 1.0 to 5 x 10^ cm“^ on a logarithmic 
scale. The white cross marks the center of the explosion. The x-axis represents the radial direction and the y-axis the direction of stellar motion (in pc). Only 
part of the computational domain is shown in the figures. Note that the panel (f) shows the supernova remnant at larger scale than in panels (a-e). 



Distance to the explosion (pc) 


Figure 7. Cross-sections taken along the direction of motion of the progen¬ 
itor through the number density field of the supernova blastwave interacting 
with the z ^ 0 region of the bow shock produced by our initially 10 Mq 
progenitor moving with velocity 40 km (Fig. [6^-1). Note that the cross- 

section from Fig. HJ) is omitted in the figure because it is very similar the 
curve (a). The respective proportion of the supernova ejecta (solid green), 
stellar wind (dotted black) and ISM material (dashed orange) are indicated 
for the curve (f). 


5 THE OLD SUPERNOVA REMNANT PHASE 

In this section we detail the interaction between supernova ejecta 
and pre-shaped circumstellar medium once the expanding front has 
passed through the bow shock. We focus on our four models of su¬ 
pernova remnants whose solutions strongly deviate from sphericity. 


5.1 Physical characteristics of the remnants 

5.1.1 Asymmetric structures... 

We show the gas density fields in the old supernova remnant 
produced by our 20 Mq progenitor moving with velocity u* = 
20 km in Fig. lAll Figs. (TO] [TT] and [12] are similar, but show¬ 
ing our 40 Mq progenitor moving with velocity u* = 20, 40 and 
u* = 70kms“^, respectively. In each of these figures, panel (a) 
corresponds to a time tysnr (our Table |4]), panel (b) shows the shock 
wave expanding into the wind cavity and panel (c) shows the rem¬ 
nant at tosnr when the simulation ends (our Table |5]). The figures 
do not show the full computational domain. In Figs. lAll tofT^ the 
overplotted solid black line is the border between the wind and ISM 
gas where the value of the passive scalar Qi{r) — 1/2, and the 
dashed black line is the discontinuity between the ejecta and the 
other materials where Q 2 (r) — 112. 

At a time tysnr, the shock wave is already asymmetric because 
its unusually dense surroundings (see our Tablejb]) strongly restrain 
it from expanding into the direction normal to the direction of mo¬ 
tion of the progenitor (panels (a) of Figs. lATl tofT^. As an example, 
the forward shock in our model OSNR2020 at a time 2400 yr has 
reached about 5.8 pc and 6 pc along the +.R and -\-z directions, 
respectively, whereas it expands to about —7.9 pc along the —z 
direction. This asymmetry of the shock wave is particularly pro¬ 
nounced in our simulation OSNR4020 whose pre-shaped circum¬ 
stellar medium contains the most massive bow shock with a mass 
of 116Mq (Fig. Ob). 

At times larger than tysnr, the shock wave freely expands into 
the undisturbed ISM in the +A and -\-z directions because it has 
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Figure 8. Density profiles (in cm“^) of our supernova remnants taken 
along the direction of motion of our initially 10 (a), 20 (b) and 40 Mq (c) 
progenitors, taken when the supernova shock wave has reached a distance 
of about 2i?(0) in the direction of motion of the progenitors, at fysnr- 


entire ly overtaken its circumstellar structure (iBrighenti & D’Ercol^ 
Il994h . see Figs. [To|[T2j)c. It progressively recovers its sphericity, 
but this takes longer in our simulations with slowly moving pro¬ 
genitors because the mass in their bow shock is larger (Fig. [Tob 
and [12]:). The penetration of the shock wave through the wake of 
the bow shock results in its chanelling into the tubular cavity of 
unshocked stellar wind (iBlandford et al.ll983h . The constant cross- 
sectional area of the cavity continues to impose large temperature 
and pressure jumps at the post-shock region of the shock wave, 
which prevents it from decelerating and which collimates the ejecta 
as a t ubular/iet-like ex tension to the spherical region of shocked 
ejecta jCox et all 199 ih . 




V-IOImVi IJlOQyf 



Figure 9. Time sequence of the evolution of the supernova remnant gen¬ 
erated by our initially 20 M© progenitor moving with velocity 20 km . 
The figures are shown at a time fysnr (a), at an intermediate time (b) and 
at the end of the simulation, at a time fosnr (c). The gas number density 
is shown with a density range from 10“^ to 5.0 cmon the logarithmic 
scale. The color scale is reversed compare to Fig. □toll] The cross marks 
the center of the explosion. Note the Napoleon’s hat structure of the unper¬ 
turbed circumstellar medium in panel (a). The solid black contours trace 
the boundary between stellar wind material and either the supernova ejecta 
or ISM gas Qi(r) = 1/2. The dotted black contour is the discontinuity 
between supernova ejecta and wind material Q 2 {r) = 1/2. The nomen¬ 
clature of the models follows our Table |5] The x-axis represents the radial 
direction and the y-axis the direction of stellar motion (in pc). Only part of 
the computational domain is shown in the figures. 
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Figure 10. As Fig. lAll with our initially 40 Mq progenitor. 


5.1.2 ...of angle-dependent physical properties 

Fig-CSplots the spatio-temporal evolution of the shock front mea¬ 
sured in our apsherical remnants along the -\-z and —z directions. 
The shock wave typically expands into the ISM at velocities of the 
order of a few hundred kms“^ whereas it propagates inside the 
trail of the bow shock with a velocity of the order of a few thou¬ 
sands km . E.g. at a time 4500 yr after the explosion, the model 
OSNR2020 has a shock wave velocity of 564 and 1578 kms“^ at 
8.73 and 13.03 pc from the center of the explosion in the direction 
along and opposite of the progenitor’s motion, respectively. The 
deceleration of the shock wave is more important for our slowly 
moving progenitors which induce the strongest anisotropy in their 
circumstellar medium (Fig. [TSh-b). whereas the ejecta velocity is 


larger after the collision with the lighter bow shocks of our fast- 
moving progenitors (Fig.fTsH and our Tab. [6]). 

Beca use the blastwa ve _ ^expand in non-uniform 

medium (iFerreira & de Jage3 l2008h . a wave created during 
the collision with the dense bow shock is reflected towards the 
center of the explosion and shocks back the unperturbed supernova 
ejecta (see Fig lAlh -c). This mechanism generates a hot region of 
ejecta which progressively Alls the entire cavity of the remnants, 
e.g. the shocked ejecta of our model OSNR2020 has n ^ 3 cm“^ 
and T ^ 10® K (Fig. lAlb c). Simultaneously, the collimated 
shock wave continues expanding downstream from the center of 
the explosion. It hits the tunnel’s side, i.e., the walls of the wind 
cavity, which post-shock density increases up to about 30 cm“® 
in model OSNR4020 and cools to less than 10® K. The sho cked 
walls produce strong optical line emission iCox et al.lll99ll. see 
also Section [ 53 ]) . 

A transmit ted wave penetrates the sho cked wind material in 
the bow shock (iBrighenti & D’Ercoi3ll994l) while a second wave 
is reflected towards the center of the tunnel. After the passage of 
the shock wave through the interface that separates the tunnel from 
the bow shock, its cross-sectional area is locally constricted and 
accelerates the flow. It is about 1500 km at a time 4200 yr, 
about 2750 km at a time 8000 yr and decelerates down to 
1450 km when the shock wave expands further in the tunnel 
at a time 12000 yr in our simulation OSNR2020 (Fig. [TSll. The 
same happens when the reflected waves collide at the center of the 
tunnel (Fig.[T3h-b). At a time tysnr almost the whole interior of the 
remnant is shocked by the reflected shock wave, and these multi¬ 
ple reflections induce a strong mixing of ejecta, wind and ISM (see 
the overlapping of the lines where Qi{r) = Q 2 {r) = 1/2). The 
Rayleigh-Taylor instabilities developing upstream from the center 
of the explosion (see, e.g. Fig.fT^) hav e an origin similar to the 
ones described in ichevalier et al.l (Il992h . but take place at the in¬ 
terface with the shocked ISM gas. Additionally, the interaction with 
the thin layer of st ellar wind niateria l elongates them up to close to 
the forward shock (iKane et alJl999t) . 

The stars end their lives as core-collapse su pernovae and the 
explo sion can produce a runaway neutron star (iLvne & Lorime j 
I 1994 I) . Assu ming a typical velo city of the compact object of about 
400kms“^ (iHobbs et al.ll2005h . one flnds that it could not be fur¬ 
ther than about 4.84, 19.8, 5.6 and 4.2 pc from the center of the 
explosion in our simulations at tosnr (Figs. IAlfc .fTob. [TTb andfTJfc). 
Consequently, we suggest that our remnants host a neutron star of 
mass Mco ~ 2 M© in the region close to the center of the explo¬ 
sion, i.e. out of the chimney-like extension of channelled ejecta. 


5.2 Remnants luminosity 

We plot the luminosities (in ergs“^) of the aspherical supernova 
remmants as a function of time (in yr) in Fig. El The bolometric 
luminosity by optically-thin radiation, L (red thin line with trian¬ 
gles), is estimated with the used cooling curve, integrating the en¬ 
ergy emitted per unit time and per unit volume over the whole rem¬ 
nant. Similarly, we plot the contribution Lism from the hot ISM 
gas (T > 10^ K, blue solid line) where radiative losses are mostly 
due to Bremsstrahlung and from the warm ISM gas (T ^ 10’^ K, 
blue thin dashed lines) where the emission are principally caused 
by cooling from Helium and metals together with O oxygen forbid¬ 
den lines emission (see Section 2.4 of Paper I) and the contribution 
Lej from the ejecta (orange dotted line). We discriminate the var¬ 
ious contributions of the luminosity L on the basis of the passive 
scalars Qi{r) and Q 2 (r) advected with the fluid. They allow us to 
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Figure 11. As Fig.[To] with space velocity Vi, = 40 km s ^. 


time-dependently trace the proportions of ISM gas, wind and su¬ 
pernova ejecta in the remnants (see section 3.4.1 and Eq. (18) of 
Paper I). The X-ray luminosity, Lx, is calculated for gas tempera¬ 
tures 10^ ^ T ^ 1.58 X 10® K, with emission coeffi cients interpo - 
lated from tables generated with the XSPEC software ( Arnaudll9^ 
with s olar m etalicity and chemica l abundances from lAsplund et al.l 
(l2009h . as in iMackev et al.l (l2015h . 

The total luminosity L is first dominated by the hot (T ^ 
10^ K) emission from the shocked gas before becoming dominated 
by emission from the warm ISM gas of temperature T ^ 10^ K at a 
time about 8000 yr (see thick solid blue line in Fig. di -d) because 
the post-shock temperature at the shock wave decreases during the 
adiabatic expansion of the blastwave. The cooling from the metals 


is stronger at T 10® K than at T ^ 10^ K (Fig.4a of Paper I), so 
L increases at larger times. 

L increases as a function of time from tysnr up to the end of 
our simulation. The reflection of the shock wave towards the center 
of the explosion produces a growing hot, dense region that is up¬ 
stream from the center of the explosion and augments the luminos¬ 
ity. The emission is infiuenced by the size of the bow shock at the 
pre-supernova phase which decreases with u* (Paper I) and governs 
the reflection of the shock wave towards the center of the explo¬ 
sion. L monotonically increases by less than one order of magni¬ 
tude over a timescale of about 10^ yr, e.g. in our model OSNR4070 
L 1.2 X 10®® ergs“^ at t 10® yr and L 9.5 x 10®® ergs“^ 
at a time 10^ yr IFig.flUi). 

Lej is smaller than Lism by a factor of a few at tysnr, e.g. just 
after the end of the shock wave-bow shock collision, L 8.0 x 
10®® and Lej 3.0 x 10®® erg in model OSNR4070 (Fig.IISi). 
It monotonically decreases with time as the shocked ejecta expands 
and its density decreases such that its emission finally becomes 
a negligible fraction of L (Fig. fT^-d). The contribution from the 
wind material is negligible compared to the bolometric luminosity. 
It is not shown in Fig. [14] since it does not infiuence L at all. 

The total X-ray emission. Lx, is calculated as the emis¬ 
sion from photons at energies 0.1-50 keV. Not surprisingly, this 
is slightly smaller than the component of Lism from the gas with 
T ^ 10^ K and follows the same trend as L (Fig. [T4k-d). The 
soft X-ray emission, i.e. from photons in the 0.5-1.0keV energy 
band, is fainter than Lx by less than an order of magnitude and 
has a similar behaviour as a function of time except for our older 
and larger supernova remnant OSNR4040 fFig.fTdbk The hard X- 
ray emission in the 2.0-5.0 keV energy band is fainter than Lx 
by abou t an order of magnitude at ty snr- It decreases as a function 
of time terighenti & D’Ercolell 19941) because the emission of very 
energetic X-ray photons ceases as the remnant expands and cools. 
Consequently, our old remnants are more likely to be observed in 
the soft energy band of X-ray emission. 


5.3 Emission maps 

In Fig. [15] we show synthetic emission maps corresponding to the 
time tosnr of the supernova remnants generated by our 20 M© 
progenitor moving with u* = 20kms“^ (panels a,d) and by 
our 40 M© progenitor moving with u* = 20 (panels b,e) and 
40kms“^ (panels c,f), respectively. Fig. im is similar for our 
40 M© progenitor moving with u* = 70kms“^. The left-hand 
side of each panel plots the Ha surface brightness (blue) whereas 
the right-hand side shows the [Olll] A 5007 surface brightness 
(green). The left-hand side of each bottom panel plots the hard X- 
ray surface brightness (red) and the right-hand side shows the soft 
X-ray surface brightness (grey). The spe ctral line emission coeffi¬ 
cients are t aken fro i n the p rescriptions bv lOsterbrock & Bochkarev! 
(Il989h and lPopit^ (Il973l) for Ha and [Qin A 5007, respectively, 
with solar oxygen abundance (Fodder^ 2003 )and imposing a cut¬ 
off temperature at 10 ® K (cf. Cox et al.lll99lh w hen oxygen be¬ 
comes further ionized ([Sutherland & Dopitalll993h . 

The region of maximum Ha surface brightness is located 
downstream from the center of the explosion. This happens be¬ 
cause the Ha emissivity, jHa oc is very faint in the region 

of hot shocked ejecta (Fig. [BJi -c). The simulation with the slowly 
moving 40 M© progenitor has an emissivity peak along the walls 
of the wind cavity (Fig. [TsT d) because effective cooling of the gas 
makes the the post-shock region cool (T > 10^ K) and dense (up 
to n ~ 50cm“®). In the other simulations, the emission origi- 
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Figure 13. Time evolution of the shock waves of our aspherical supernova remnant models. The figure plots the shock wave position (thin lines, in pc) and 
velocity (thick lines, in kms“^), measured along the direction of motion of the progenitor, both upstream (dashed red lines, -\-z direction) and downstream 
(solid blue lines, —z direction) from the center of the explosion, respectively. The figure represents the expansion from fysnr up to the end of the simulation at 
fosnr (in yr) of the remnants generated by our initially 20 Mq progenitor moving with velocity 20 km (a), our initially 40 Mq progenitor moving with 
velocity 20 (b), 40 (c) and 70 km (d), respectively. 


nates from the outer layers of the bow shocks because the walls of 
their less massive bow shocks allow a faster and deeper penetra¬ 
tion of the shock wave into the shocked wind material (Fig. 

Our predicted maximum Ka emission is above the diffuse emis¬ 
sion sensitivity limit of the SuperCOSMOS H-alpha Survey (SHS 


IParker et 'Zl l2005h of 1.1-2.8 X 10“ ^^ergs ^ cm ^ arcsec 
e.g. our model OSNR4040 has a maximum emission larger than 
9 X 10~^® ergs“^ cm“^ arcsec“^ (Fig. [Bh), and could be com¬ 
pared with data from these surveys. 

The [Olll] maximum surface brightness of our models origi- 
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Figure 14. Lightcurves of our old apherical supernova remnants. We show the luminosities (in ergs“^) during the time interval [tysnr,fosnr] (in yr) of 
our aspherical supernova remnants generated by our initially 20 M© progenitor moving with velocity 20 km (a), our initially 40 M© progenitor moving 
with velocity 20 (b), 40 (c) and 70 km (d), respectively. The hgures distinguish between the bolometric (total) luminosity from the gas (red thin line with 
triangles), the contribution from the hot ISM gas (T > 10^ K, blue solid line), from the warm ISM gas (T ^ 10^ K, blue thin dashed lines) and the contribution 
from the ejecta (orange dotted line). The black lines are the X-ray luminosity calculated at photons energies E > 0.1 keV (T > 1.2 x 10® K, dashed line), 
in the soft energy band 0.5-1.0keV (5.8 x 10® < T < 1.2 x 10^ K, crossed line) and hard energy band 2.0-5.0keV (2.3 x 10^ < T < 5.8 x 10^ K, 
squared line), respectively. 


nates from the dense (n ~ 2—10 cm ^) and warm (T ^ 10® K) 
post-shock region behind the shock wave. It is located at the walls 
of the cavity and produces a ringed/tubular structure (Fig. [Bt-c 
andfT^). It is generally not coincident with the projected Ha emis¬ 
sion because the [Olll] emissivity has a different dependence on 
the temperature, i.e. j[oiii] cx: However, the simula¬ 

tion OSNR4020 with a heavy bow shock have its [Olll] maximum 
surface brightness all along the walls of the wind tunnel and it is 
coincident with the region of maximum Ha emission, i.e. the be¬ 
haviour of the emissivities with respect to the large compression of 
the gas in the walls (oc v?) overwhelms that of the gas temperature 
(Fig.[l5]3). 

The X-ray emission originates from the hot gas with 10^ ^ 
T ^ 10® K, i.e. from near the shock wave expanding into the un¬ 
perturbed ISM. The maximum surface brightness comes from the 
hot region of doubly shocked gas upstream from the center of the 
explosion, and from the post-shock region at the channelled shock 
wave (Figs. [T5h-f.[Tbb). The hard X-ray surface brightness is sev¬ 
eral orders of magnitude fainter than the projected soft X-ray emis¬ 
sion, because the gas is not hot enough at tosnr (Fig.O. The emit¬ 
ting region in soft X-rays is peaked in the post-shock region at the 
shock wave whereas the hard X-ray come from a broader region 
of shocked gas which outer border is the shock wave (Fig. 

Note also the anti-correlation between the surface brightness in X- 


ray (hot region from near the forward shock) and in [Olll] (colder 
and denser walls) in our remnants generated by a fast-moving 
(u* ^ 40kms“^) progenitor (Fig. [I5]3,e;c,f). 


6 DISCUSSION 

We here discuss our results in the light of precedent studies and 
pronounce on the best manner to observe aspherical supernova rem¬ 
nants generated by massive Galactic runaway stars. Finally, we ex¬ 
amine our models in the context of recent observations. 

6.1 Comparison with previous works 

We tested that the code PLUTO (iMignone et al.l l2007l . l2012l) re¬ 
produced the one-dimensional models of core-collapse super- 
nov ae interacting with their surroundings in I Whalen et al.l (l2008l) 
and Ivan Veelen et aP (l2009l) using a uniform, spherically symmet¬ 
ric grid_Ou£jmmeriuaI method (Paper I) is different from that 
in I Whalen et al.l (l2008l) because (i) they utilise a finite-difference 
scheme coupled to a network of chemical reactions following 
the non-equilibrium rates of the species composing the gas and 
(ii) their algorithm inc ludes artificial viscosity (ZEUS code, see 
IStone & Norman|[l992h . 
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Figure 15. Synthetic emission maps of our models OSNR2020 at a time 12100 yr (a,d), OSNR4020 at time 49500 yr (b,e) and OSNR4040 at a time 
14000 yr (c,f), plotted on the linear scale in units of ergs“^ cm“^ arcsec“^. The left-hand part of the figures show the Ha surface brightness (blue) and 
the right-hand part the [Oiii] surface brightness (green), respectively. The black crosses mark the center of the explosion. The panels d-f show the projected 
X-ray emission maps in the hard (2.0-5.0keV corresponding to 2.3 x 10^ < T < 5.8 x 10^ K, left-hand part) and soft (0.5-1.0keV corresponding to 
5.8 X 10® < T < 1.2 X 10^ K, right-hand part) X-ray bands for the same models. The a:-axis represents the radial direction and the y-axis the direction of 
stellar motion (in pc). Only part of the computational domain is shown in the figures. 


W e ran tests with their cooling curve (iMacDonald & Bail^ 
Il98lh and with a cooling curve for collisional ionization equi¬ 
librium medium (see details in section 2.4 of Paper I). We find 
no notable differences, mostly because they are similar in the 
high temperature regime (T ^ 10^ K) that is relevant for the 
supernova-wind interaction (Fig.|5]). We extend this method to two- 
dimensional, cylindrically symmetric tests of a supernova shock 
wave expanding into a homogeneous ISM to ensure that the 
sphericity of the shock wave is conserved throughout its expan¬ 
sion. We notice that the solution behaves well with respect to the 
symmetry axis Oz. 


Models of an off-centered explosi on in a wind-driven cavity 
are available in iRozvczka et al.l (Il993h . Their model produces a 
parsec-scale jet-like feature as do our aspherical models (Figs. lAlb 
to m) but they neglect the progenitor’s stellar evolution, as¬ 
sume a different microphysics, and correspond to a totally differ¬ 
ent point of the parameter space (n = 10® cm“®). Our descrip¬ 
tion of the supernova shock wave interacting with its pre-shaped 
bow shock is consistent with the wo rks tailored to the Kep ler’s 
sup ernova remnant; s ee sec t ion 5.2 of Borkowski et al.l (Il992h but 
also I Velazquez et ^(l2006h : lToledo-Rov et al.l(l2014l) . 


The tunnel of unshocked wind that channels the shock wave. 
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Figure 16. As Fig. El with an initially 40 M© progenitor moving with space velocity r;* = 70 km s ^. The Ha and [Oiii] A 5007 surface brightnesses are 
plotted in the left-hand panel, the projected hard and soft X-ray emissivity in the right-hand panel. 


e.g. in our model OSNR4020 (Fig. lAlh ). is morphologically and 
structurally consistent with the model of the red supergiant pro¬ 
genitor of the Crab nebula in ICox et^ (Il99lh . assuming a diluted 
ambient medium and a larger space velocity (n = 0.25 cm“^, 
u* = 69.5 km s“^). Notice that the simulations with fast-moving 
stars or post-main-sequence high-mass-loss events have a tunnel 
with clumpy walls (Fig. lAlh and \T2k ) which do not p revent the 
tunnel ing of the shock, in contrast to the suggestion by ICox et al.] 
(Il99lh . 

The growth and overall shape of our aspherk al remnants 
are in accordance with iBrighenti & D’Ercold (1 19941) . Our model 
OSNR2020 is morphologically consistent with their model 1 (u* = 
17kms“^, n = 1 cm“^). They assumed a comparable mass loss 
(M = 10“^Moyr“^) but a twice larger wind velocity (uw = 
20kms“^) during the red supergiant phase which lastes about 
7 X 10^ yr. Identical remarks arise comparing their model 3 and 
our simulation OSNR4070 (Fig. [T2^-c). Because our models in¬ 
clude thermal conduction, the region of shocked wind in the bow 
shock from the main-sequence phase is larger and the tunnel in our 
simulation OSNR4020 narrower than in their model 4 and allows 
a more efficient channeling of the shock wave (Fig.[T0k-c). More 
details about the effects of heat conduction in our remnants is given 
in Appendix IaI 

As a conclusion, for overlapping parameters our results agree 
well with previous models of supernova remnants produced by run¬ 
away progenitors. We extend the parameter space with a repre¬ 
sentative sample of models tailored to the Galactic plane, whose 
progenitor’s wind properties are ta ken from self-con sistently pre¬ 
calculated stellar evolution models (iBrott et al.ll201lh . 


6.2 Comparison with observations 

6.2.1 General remarks and comparison approach 

Comparing our simulations with observations is not a straight¬ 
forward task. Even though this paper explores a representative 
sample of Galactic, massive, runaway stars, our remnants can 
only develop asymmetries when the isotropic shock wave inter¬ 
acts in reality with a dense bow shock, whereas other mecha¬ 
nisms can also induce asymmetrie s. They can originat e from an 
intrinsically anisotropic explosion (iBlondin et al.lll99^ . the rota¬ 


tion of the progenitor (jL anger et al.l 19991) . the magnetization of the 
ISM (IRozvczI^ & TenOTio Taglell 19951) or the presence of a neigh- 
bouring circumstellar structure (iFerreira & de Jagei]l2008l) . Conse¬ 
quently, we hereby simply attempt to establish a qualitative dis¬ 
cussion between the density fields in our models and observations 
available in the literature. 

Su pernova reni nants are often observed in radio wave¬ 
lengths (lGreenll200^ . e.g. at about 325, 843 or 1400 Mhz. This 
emission arises from bremsstrahlung (proportional to n^, where Ue 
is electron number density ), synchrotron radiation (proportional 
to UeB, where B is the field strength), and maser emission (tracing 
very dense ISM regions). Our discussion is limited to a simple com¬ 
parison between such published measures and the density fields 
in our hydrodynamical simulations. Additionally, because Galactic 
SNRs usually have poorly constrained distances (and hence sizes 
and masses), we discuss morphological similarities even though it 
may not be clear whether or not our model parameters are appro¬ 
priate to a given observed remnant. 

6.2.2 Observations from maser emission such as in 3C391: 

explosions in a dense medium, runaway progenitors or 
both ? 

IBrighenti & D’Ercc^ (1 19941) and lEldridge et al.l (1201 ll) already 
justified the relevance of studying runaway O stars in the un¬ 
derstanding of supernova remnants and gamma-rays bursts. Par¬ 
ticularly, they underline the difficulty of interpreting the shape 
of incomplete and/or inhomogeneous arc-like supernova remnants 
because the overdensities upstream from the center of the ex- 
plosion can slIso ar ise from pre-existing dense regions (see, e.g. 
lOrlando et al.ll200i) . The presence of OH maser emission orig¬ 
inating fr om the shock front can discriminate be tween these 
scenarios (lErail et al.l 1 19961 : lYusef-Zadeh et al.l l2003l) . and led to 
the classification of about 20 Galactic arced remnants such as 
G31.9-1-0.0 or G189.1+3.0 as running into a dense cloud. 

Examining the necessary conditions to produce OH 1720 Mhz 
(T 50-125 K, rimedi iim ^ 10^ cm OH column density 
Nou ^ 10^® cm“^, see iLockett et al.l[l9^ one finds that such 
high densities are not reached in our models. However, accord¬ 
ing to the Rankine-Hugoniot jump relations it may be the case 
in the post-shock region if the shock propagates in a medium 
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of pre-shock density nmedium/4 w 2.5 x 10"^ cm“® such as a 
starless dense core, a molecular cl oud or a contracting ISM fil¬ 
ament (iKaufman & Kaufm"^ l2009h . Assuming a temperature of 
T 50K, the medium sound speed is Cs ~ 0.65 kmso any 
star moving with velocity u* ~ 5cs ~ 3.25 km would have 
an hypersonic motion and produce a bow shock. Consequently, if 
OH maser emission is not likely to be detected from bow shocks in 
the field, it naturally arises from the surroundings of runaway stars 
in a dense medium. Moreover, if the star explodes in such an en¬ 
vironment, additional constraints may be necessary in order to dis¬ 
tinguish between maser emission produced because of the runaway 
nature of the progenitor, emission originating from an expanding 
shock wave, or both. 

Particularly, 3C391 was originally beli eved to be the archety¬ 
pal r emnant from a moving progenitor (iBrighenti & D’Ercola 
I 1994 I) but is now known to be associated with OH maser emis¬ 
sion (iFrail et al.l Il996h together with molecular line emission. 
The Infrared Spatial Observator y (ISO) revea l ed, am ong other, 
CO, HCO^ line emis sion (see Reach & Rhol 1 19991 . and refer¬ 
ences therein) whereas iNeufeld et al. ( 2014 ) reported recent in¬ 
frared analysis of the H 2 and H 2 O line emission. However, as men¬ 
tioned above, nothing indicates whether this emission originates 
from the shock wave colliding with a wind bubble produced by 
the progenitor or from the shock wave running into its dense sur¬ 
rounding medium. Modelling the explosion of a supernova from a 
runaway progenitor with rimedium ^ 1 cm“^ is conceivable but 
beyond the scope of the present investigation and left for a follow¬ 
up study. Considering the past evolution of supernova progenitors 
in the modelling of their remnants will help to discriminate between 
these two situations. Further investigations, e.g. tailored to 3C391, 
may provide more severe constraints on its circumstellar medium 
at the pre-supernova phase and help to understand its formation 


as the field is strong enough (iRozvczka & Tenorio-Taglelll995h and 
produ ce X-ray and/or radio synchrotron emission f r om the opposed 
arcs (IVelazquez et^ 120041 : IPetruk et ^ l2009l : ISchneiter et al.l 

l2010h . 


6.2.4 The jet/tubular-like extension shaped by the motion of the 
progenitor 

Our Galactic, slowly moving, initially 20 M© progenitor produces 
a supernova remnant whose outer region strongly emits in [O III]. 
The remnant generated by our slowly moving, initially 40 M© 
progenitor has an [Olll] jet-like feature that has a Ha counter¬ 
part (Fig. (TS]) generated by ejecta channelled with velocity about 
1000 km into the wind tunnel of the bow shock (Fig. \T3h). 
Those tubular/jet-like features (Figs. [T5l[T^ are re miniscent of the 
chimn ey discovere d in [Olll] in the Crab nebula JBlandford et al.l 
Il983h modeled in ICox et ah (Il99lh . This morphological resem¬ 
blance mostly arises because of the similar stellar evolution history 
and space velocity u* 30 km in both simulations. 

Nevertheless, further investigations and detailed comparison 
would have to take into account the youth of the Crab nebula 
1000 yr) and its location in the high latitude, low-density ISM 
of the Galaxy. In our simulations with larger u*, those jet-like 
extensions become an [O III] tubular structure that is thinner and 
closer to the throttling separating the surroundings of the center of 
the explosion from the trail of the bow shock. Supernovae explod¬ 
ing in a wind cavity could hence form tunnels or barrel-like shapes, 
however, alternative quite convincing demonstration, e.g. based 
on asymmetric explosions hav e already been proposed (WB49, 
iGonzalez-Casanova et al.ll2014l) . 


6.2.3 Bilateral supernova remnants 

Our simulations PSN4020 and PSN4040 produce approximately 
cylindrical wind-blown cavities with dense walls. The interaction 
of the supernova with this circumstellar medium, in our models 
OSNR4020 and OSNR4040 (Figs. [T0| and [TT]) produces the con¬ 
ditions required for a bilateral supernova remnant, i.e. a cylindrical 
cavity with dense, shocked gas on the boundaries. A previous study 
identified the material composing the bilateral struc tures of the 
remn ant G296.5-1-10.0 as shocked pre-supernova wind dManchesteii 
Il987h and the prese nce of a neu t ron st ar in between the opposed 
arcs is reported in IZavlin et^ (Il998h . It constrains the stellar 
remnant to be a 1.4 M© neutron star, so its pro genitor had initial 
mass below about 25 M© (IWooslev et al.ll2002h . It is similar to our 
slowly moving model of an initially 40 M© at time 11400 yr after 
the explosion, when the shock wave hits the very dense walls of the 
wind tunnel (Fig.fTOb). Our understanding of the bilateral charac¬ 
ter of this remnant as a result of its progenitor’s supersonic motion 
only is valid if one considers a lower initial mass progenitor mov¬ 
ing in a medium dense enough to form heavy cavity walls resistant 
to the shock wave. No te that this result is consistent with fig. 6D 
of lOrlando et al.l (l2007h . 

However, alternative explanations have been proposed for the 
fomation of bilateral remnants like G327.6-1-14.6 or G3.8-0.3. A 
strong axisymmetric background ISM magnetic field (not included 
in our simulations) has also been suggested to be responsi ble for the 
bilate ral character of some Galactic supernova remnants (iGaenslej 
Il998h . It would indeed make their shape more elongated as long 


6.2.5 An alternative explanation for the Cygnus loop nebula ? 

G074.0-8.5, also called the Cygnus loop, is a supernova remnant for 
which X-ray observations, e.g. with ROSAT (lAschenbach & Leah^ 
Il999h reveal a characteristic overall drop-like morphology. An 
early interpretation of this s hape is a champagne blow out from the 
edge of a molecular cloud (iTenorio-Tagle et al]ll985h . Recent X- 
ray data favours a model with an explosion into a pre-shaped cav- 
ity and whose ejec ta have recently impacted the imperfect walls. 
lUchida et al.l (l2009l) supports this model and derives the remnant’s 
age to be about 10000 yr, its radius to be 12-17 pc and its initial 
progenitor’s mass to be 12-15 M©, which is consistent with our 
model OSNR2020 (Fig.[T5]i). 

Moreover, some southern optical filame nts have a measured 
expansion velocity of a few hundred kms“^ (iMedina et al.ll2014l) 
as in our simulation OSNR2020 (Fig.[T3h). In our picture, the south 
blowout region of the remnant is the low-density wake behind the 
pr ogenitor in which the gas is the hottest (> 10® K), as noted 
in lAschenbach & LeahvI (Il999h . Nevertheless, our models predict 
that the soft and hard X-rays emission should be limb-brightened, 
whereas the observations have a filled morp hology. This may be 
explained by the presence of a neutron star jKatsuda et al.l 120121) 
that is not taken into account in our models. Future simulations, 
investigating the plerionic nature of the Cygnus loop nebula will 
allow us to assess this interpretation. 
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6.2.6 Distinguishing remnants from moving progenitors and 
remnants with ISM-induced anisotropy. 

Distinguishing whether a one-sided supernova remnant is produced 
thanks to the motion of its progenitor or because of its surround¬ 
ings’ inhomogeneity is difficult since both situations can happen to¬ 
gether. As discussed above, the OH emission can trace the presence 
of a molecular cloud, but is not sufficient to exclude the presence of 
an hypothetical bow shock. A couple of additional comments can 
be given. 

(i) In the situation of a runaway supernova progenitor, the cavity 
that channels the subsequent shock wave left behind its bow shock 
may produce a more collimated remnant than in the case of an in¬ 
teraction with a dense cloud. This phenomenon could be enhanced 
by the presence of a background ISM magnetic field which elon¬ 
gates the pre-shaped tunnel, reducing the instabilities affecting its 
walls and consequently favoring the channeling. 

(ii) If a runaway massive star sheds (an) heavy envelope(s) 
throughout its evolution, e.g. via a luminous blue variable event or a 
Wolf-Rayet phase, the passage of the shock wave through these ex¬ 
pelled shell(s) will fragment them. Their temperature will increase 
by shock heating before cooling down, inducing an X-ray rebright- 
ning of the formed foculli. Such a mixed-morphology supernova 
remnant generated by a blown bow shock would have an incom¬ 
plete arc-like radio envelope around a central region containing X- 
ray-emitting clumps. 

(iii) The situation becomes more complicated when the explo¬ 
sion happens at the interface between two different phases of the 
ISM. If a runaway progenitor leaves a molecular cloud, the tubular 
extension produced by the channeled shock wave will be directed 
towards the center of the cloud. If the progenitor enters the cloud, 
the imprint of the last mass-loss events of the moving star onto 
the border of the cloud will produce a characteristic hole. A fa¬ 
mous (nonetheless extragalactic) example of such a phenomenon is 
the structure that o verhangs the rem nant of SN1987a in the Large 
Magellanic Cloud JWang et al.lll993h . 


7 CONCLUSION 

In this paper, we present a grid of hydrodynamical models of asym¬ 
metric supernova remnants generated by a representative sample 
of Galactic runaway massive stars whose circumstellar medium 
d uring the ma i n-sequ ence and red supergiant phases is studied 
in iMever et al.l (l20l4) . We compute the bow shocks generated by 
our progenitors from near the pre-supernova phase and model the 
collision between the supernova shock waves and the circumstellar 
medium which result in the generation of supernova remnants. The 
progenitors’ initial masses range from 10 to 40 M© and they move 
with space velocities ranging from 20 to 70kms“^. Our models 
include both optically-thin cooling and photoheating of the gas. 
Electronic thermal conduction is included in the calculations of the 
circumstellar medium and in the simulations of the supernova rem¬ 
nants. 

We stress that the stellar motion of a core-collapse supernova 
progenitor can be responsible for the deviations from sphericity of 
its subsequent remnant. We show that the bow shocks trapping at 
least 1.5 M© of ISM gas are likely to generate aspherical super¬ 
nova rem nants. They correspond to high mass and/or slow ly mov¬ 
ing stars terighenti & D’Ercolelll994riMever et al.ll20l3l . At the 
ISM number density that we consider, they are produced either 
by our initially 20 M© star moving with space velocity of about 


20kms“^ or by our initially 40 M© runaway star. These mass- 
accumulating bow shocks generate a dense bulge of shocked ISM 
gas upstream from the direction of motion of the star whereas a 
cavity of low-density wind material forms in the opposite direction 
and extends as a tunnel of unshocked wind material into the trail of 
the bow shock. 

After the supernova explosion, the shock wave expansion is 
strongly infiuenced by the anisotropy of its circumstellar medium. 
It collides with the overdense part of the bow shock whereas it ex¬ 
pands freely at velocities of the order of 1000 km in the oppo¬ 
site direction, channelled by the pre-shaped tu nnel of unshocke d 
wind material (cf. observations of RCW86 in I Vink et~^ Il997h . 
The mass of shocked ISM trapped in the bow shock decelerates 
the shock wave, which continues to penetrate the unperturbed ISM 
after the collision with velocity of the order of 100 km . 

As the shock wave evolves in a non-uniform medium, 
it is partly refiected towar ds the center of the explo¬ 
sion (lEerreira & de Jage3 l2008h after the collision with a dense 
bow shock. It induces mixing of supernova ejecta, stellar wind and 
ISM gas that is particularly important for fast-moving progenitors. 
This refiected wave shocks the zone around the center of the 
explosion, reheating the gas, which subsequently cools below 
10^ K because of the adiabatic expansion of the blastwave. Its 
luminosity increases, dominated by thermal Bremsstrahlung and 
soft X-ray emission originating from the shocked ISM that is 
upstream from the center of the explosion. The emission from the 
ejecta or from the progenitor’s wind material does not contribute 
significantly to the remnants’ total luminosity once the bow shock 
is overtaken by the supernova shock wave. 

Our Galactic aspherical supernova remnants have an [Olll] 
A 5007 surface brightness larger than their projected Ha emission, 
i.e. the [Olll] is more appropriate line to search for Galactic su¬ 
pernova remnants. Their [Olll] surface brightness is maximum in 
the post-shock region of the shock wave. It is concentrated along 
the walls of the tunnel of wind material. The region of maximum 
Ha emission is downstream from the direction of motion of the 
progenitor. It originates from the outer part of shocked ISM ma¬ 
terial in the trail of the progenitor’s bow shock. In the case of our 
slowly moving initially 40 M© progenitor, it mainly comes from 
the region where the shock wave interacts with the walls of the tun¬ 
nel, i.e. the ejecta forms an [O III] tubular-like feature that has an 
Ha counterpart. Moreover, we find that our remnants are brighter 
in soft X-ray emission originating from near the shock wave than 
in hard X-ray emission coming from the post-shock region at the 
shock wave. 

Supernova remnants generated by runaway progenitors natu¬ 
rally show structures highly deviating from sphericity. Particularly, 
our models of remnants generated by high-mass, slowly moving 
progenitors have morphologies consistent, e.g. with the bilateral 
character of observed barrel-like Galactic supernova remnants such 
as G296.5-I-10.0 or with the morphology of the Cygnus loop neb¬ 
ula. However, other mechanisms are at work in the shaping of su¬ 
pernova remnants. Eorthcoming work will investigate the effects of 
an ISM magnetic field on the evolution of our remnants, in order 
to quantitatively appreciate its consequences on the remnants’ dy¬ 
namics and emission signatures. 
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Al The effects of thermal conduction on the pre-shaped 
cirumstellar medium 

In the situation of our slowly-moving stars (u* = 20 km s“^), we have an 
off-centered explosion inside the wi nd-bubble generated during the main- 
sequence phase of the moving star (iBrighenti & D’Ercolgll994l) . The ef¬ 
fects of thermal conduction mainly arise from the structural differences it 
induces in the bow shock, which consist in the circumstellar medium of the 
star at the pre-supernova phase (see our model PSN4020 in Eig. lAlh k The 
high effective temperatures of our m ain-sequence, massiv e stars produce a 
heat flux which timescale oc (lOrlando et al.ll2005L Raper I) is much 

faster than both the dynamical advection timescale of the stellar wind and 
ISM gas in the bow shock and cooling by optically-thin radiative processes 
timescale. It transports large amount of the gas internal energy from the re¬ 
verse shock of the bow shock towards the contact discontinuity that slightly 
enlarges the bow shock in the direction of motion of the star. Moreover, it 
damps the instabilities developing at the reverse shock, i.e. the walls of the 
tunnel, and reorganises its internal structure (see discussion in section 3.3 
of Raper I). 

In the case of our fast-moving stars (u* ^ 40kms“^), we have an 
off-centered explosion outside the wind-bubble generated during the main- 
sequence phase of the moving star and the supernova explosion conse¬ 
quently takes place in a pre-shaped circumstellar medium where any in- 
for mation relative to the rnain-se quence wind bubble is located behind the 
star (IBrighenti & D’Ercolel 19941) . The effects of thermal conduction mainly 
arise from the instabilities that develop in the bow shock produced during 
the post-main sequence phase of the stellar evolution. Since the seeds and 
the growth of the Kelvin-Helmholtz and Rayleigh-Taylor instabilities partly 
depend on density distribution at the end of the main-sequence phase, the 
effects of thermal conduction consist in changing the shape of the instabili¬ 
ties at the apex of these bow shocks (see our model RSN4070 in Fig. lAlb k 
Note that only the morphology of the eddies at a given time differ, the de- 
velopi ng instability remains of the same kind, i.e . a thin-shell-related insta¬ 
bility (IVishniaclll9'94lBlondin & KoerweJl998h . 

A2 The effects of thermal conduction on the old supernova 
remnants 

In the situation of our slowly-moving supernova progenitors, the remnant 
is affected by thermal conduction in the sense that the channeled supernova 
shock wave interacting with the walls of the cavity encounters a more irreg¬ 
ular and clumpy medium (Eig.|A2t). It changes neither the shape of the out¬ 
flow along the direction of the progenitor nor the global remnant’ morphol¬ 
ogy, however, it may generate ring-like emission artefacts, e.g. by projection 
effect. Note that heat conduction i s in its turn sensitive to the ISM magneti¬ 
sation which makes it anisotropic dBalsara et al.l200^ as it has been shown 
in several studies de voted to old supernova remnants d Velazquez et al .12004 
iBalsara A alJ2008h . 

In the case of our fast-moving supernova progenitors, the absence of 
thermal transfers (i) slightly modify, as discussed above, the shape of the ed¬ 
dies constituting the pre-supernova circumstellar medium distribution and 
(ii) allows the ejecta interacting with the bow shock to cool independently of 
the rest of the gas during the radiative phase of the remnant. Consequently, 
the Angers of the developing Rayleigh-Taylor instabilities become denser 
(Fig. EH)). 


APPENDIX A: THE EFFECTS OF THERMAL 
CONDUCTION 


In this Appendix we discuss the effects of thermal conduction on the shape 
and internal structures of our circumstellar nebulae, before and after the 
supernova explosion. 
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Figure Al. Effects of heat conduction on the pre-supernova circumstellar 
medium of our initially 40 Mq progenitor moving with velocity 20 km s“ ^ 
(a) and 70kms“^ (b). The figure shows the gas number density field of 
simulations carried out without (left-hand part of the panels) and with heat 
conduction (right-hand part of the panels) with a density range from about 
10“^ to 5.0 cmon the logarithmic scale, at a time tosnr- The white 
cross marks the position of the runaway star. 


Figure A2. Effects of heat conduction on the old supernova remnant gen¬ 
erated by our initially 40 Mq progenitor moving with velocity 20 km 
(a) and 70kms“^ (b). The figure shows the gas number density field of 
simulations carried out without (left-hand part of the panels) and with heat 
conduction (right-hand part of the panels) with a density range from about 
10“^ to 5.0 cm“^ on the logarithmic scale, at a time tosnr. The white 
cross marks the center of the explosion. 
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